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Abstract 

We present an improved measurement of the Rossiter-McLaughhn effect for the ex- 
oplanetary system HD 17156, based on radial- velocity data gathered with the Subaru 
8.2m telescope throughout the planetary transit of UT 2008 November 7. The data 
allow for a precise and independent determination of the projected spin-orbit an- 
gle of this system: A = 10.0° ±5.1°. This result supersedes the previous claim of 
A = 62° ±25° by Narita et al., which was based on lower-precision data with poor 
statistics. Thus the stellar spin and planetary orbital axes of the HD 17156 system 
are likely to be well-aligned, despite the planet's large orbital eccentricity suggesting 
a history of strong dynamical interactions. 
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Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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1. Introduction 

The discoveries of more than 340 extrasolar planets have revealed the remarkable diver- 
sity of planetary systems. For example, many Jovian planets lie far inside the snow line, and 
many of their orbits are more eccentric than the orbits in the Solar system. This diversity of 
planetary systems has prompted theorists to consider new modes of planet formation and mi- 
gration, including Type II migration models (e.g., Lin & Papaloizou 1985; Lin et al. 1996; Ida & 
Lin 2004), planet-planet scattering models (e.g., Rasio & Ford 1996; Marzari & Weidenschilling 
2002; Nagasawa et al. 2008; Chatterjee et al. 2008), and Kozai cycles (with tidal friction) for 
planets in binary systems (e.g., Wu & Murray 2003; Takeda & Rasio 2005; Fabrycky & Tremaine 
2007; Wu et al. 2007). In particular, planet-planet scattering and the Kozai effect may help to 
explain the eccentric orbits of close-in planets, since Type II migration alone does not seem to 
produce large orbital eccentricities. 

Observations of transiting planets may lead to better understanding of migration mech- 
anisms through the Rossiter-McLaughlin effect (hereafter the RM effect: Rossiter 1924, 
McLaughlin 1924). The RM effect is an apparent radial- velocity anomaly that is observed 
during planetary transits (see Ohta et al. 2005; Gimenez 2006; Gaudi & Winn 2007 for theo- 
retical discussions). By measuring and modeling this anomaly, one can learn the sky-projected 
angle between the stellar spin axis and the planetary orbital axis, denoted by A. The angle 
between these two axes has been recognized as an useful diagnostic to investigate planet migra- 
tion histories. The basic premises are that planets are formed with their orbital axes parallel to 
the stellar rotation axis, that Type II migration maintains a small spin-orbit alignment angle, 
while planet-planet scattering and the Kozai effect may produce large spin-orbit misalignments 
(e.g., Nagasawa et al. 2008; Chatterjee et al. 2008), and that tidal effects do not subsequently 
alter the spin-orbit angle (Barker & Ogilvie 2009). To the extent that these premises are valid, 
observations of the RM effect in transiting planets allow us to test planet migration models. 
One may compare the distribution of observed spin-orbit alignment angles with theoretical 
predictions based on the scenarios of planet-planet scattering or Kozai cycles with tidal friction 
(Fabrycky & Winn 2009). 

Measurement of the RM effect is especially interesting for transiting planets in eccentric 
orbits (hereafter "transiting eccentric planets" or TEPs), because the migration mechanisms 
that produce eccentric orbits may also produce large spin-orbit misalignments. The number of 
known TEPs is still small. Observations of the RM effect have been reported for a few cases, 
including HAT-P-2b (Winn et al. 2007b; Loeillet et al. 2008), X0-3b (Hebrard et al. 2008; Winn 
et al. 2009b), and the subject of this paper, HD 17156b. 

HD 17156b was discovered in a Doppler survey (Fischer et al. 2007; F07) and later 
found to transit its parent star (Barbieri et al. 2007). This planetary system stands out among 
the transiting planets for its large orbital eccentricity (e = 0.68) and relatively long orbital 



2 



period (21 days). The host star is very bright {V = 8.2), facihtating the detection of the RM 
effect and many other follow-up observations. For the HD 17156 planetary system, previously 
Narita et al. (2008; N08) reported a possible spin-orbit misalignment A = 62° ± 25°, based 
on radial- velocity data obtained at the Okayama Astrophysical Observatory (OAO). A more 
decisive conclusion, however, required further precise observations with better statistics. Indeed 
Cochran et al. (2008; C08) and Barbieri et al. (2008; B08) obtained new radial-velocity data 
(C08: HET and HJST datasets, B08: TNG dataset), and claimed a close alignment, A = 
9.4° ± 9.3° (C08) and A = -4.8° ± 5.3° (B08)^ 

We report here on our new observations with the Subaru 8.2m telescope, covering a full 
transit of HD 17156b and providing the larger and more precise radial- velocity dataset for the 
system around the transit phase. Section 2 summarizes our Subaru spectroscopic observations 
and radial- velocity dataset, and section 3 describes the analysis procedure for the RM effect. 
Section 4 presents and discusses the results for the derived system parameters. Section 5 
summarizes the main findings of the present paper. 

2. Subaru Spectroscopy and Radial Velocities 

We observed HD 17156 over 6.5 hr covering a full transit of HD 17156b with the High 
Dispersion Spectrograph (HDS: Noguchi et al. 2002) on the Subaru 8.2m telescope on UT 2008 
November 7. We employed the standard I2a set-up of the HDS, covering the wavelength range 
4940 A < A < 6180 A and using the Iodine gas absorption cell for precise radial velocity mea- 
surements. The slit width of 0''4 yielded a spectral resolution of ~90000, and the seeing on 
that night was between O'.'A and 0''7. The exposure time was 3-5 minutes, yielding a typical 
signal-to-noise ratio (SNR) of approximately 180 per pixel. We processed the observed frames 
with standard IRAF^ procedures and extracted one-dimensional spectra. We computed rela- 
tive radial velocities following the algorithm of Butler et al. (1996) and Sato et al. (2002), as 
described in Narita et al. (2007). We estimated the internal error of each radial velocity as the 
scatter in the radial-velocity solutions among the ~4 A segments of the spectrum. The typical 
internal error is 3 m s~^. The radial velocities and their internal errors are plotted in Figure 1 
and summarized in Table 1. 

3. Analysis Procedure 

We empirically model the RM effect of HD 17156 following the procedure of 
Winn et al. (2005). Specifically, we simulate spectra that are affected by the RM effect, and 

^ Note that BOS used /3(= —A) defined in Gimenez (2006) for the projected spin-orbit ahgnment angle. 

^ The Image Reduction and Analysis Facility (IRAF) is distributed by the U.S. National Optical Astronomy 

Observatories, which arc operated by the Association of Universities for Research in Astronomy, Inc., under 

cooperative agreement with the National Science Foundation. 
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then we "measure" the anomalous radial velocity of the simulated spectrum using the same 
algorithm that we use on the actual data. We start with the NSO solar spectrum (Kurucz et al. 
1984) and deconvolve it to remove the rotational broadening of the Sun. We neglect the differ- 
ential rotation of the Sun and adopt 1.85 km s~^ (the rotational velocity at the equator) for the 
deconvolution kernel. We then apply a rotational broadening kernel with Vsinlg = 2.6 km s~^ 
as appropriate to mimic the disk-integrated spectrum of HD 17156 (F07). In this step, for 
simplicity, we neglect differential rotation, microturbulence and macroturbulence, convection 
cells, and possible starspots and other active regions. Finally, to simulate the RM effect, we 
subtract a scaled and velocity-shifted copy of the original unbroadened spectrum, representing 
the portion of the stellar disk hidden by the planet. We produce many simulated spectra using 
different values of the scaling factor / and the velocity shift Vp, and compute the apparent 
radial velocity of each spectrum. We thereby determine an empirical relation between the size 
of the planet and its position on the stellar disk, and the apparent radial velocity of the star. 
We fit a polynomial function to this relation, and find 



As has been found previously (e.g., Winn et al. 2005), this formula differs from the analytic 
formula presented by Ohta et al. (2005). The two formulas are nearly consistent, however, if 
we rescale the value of l^sinJs in the analytic OTS formula by a factor of 1.2. In this sense 
the OTS formula may overestimate the projected stellar rotational velocity by ~20%. Further 
discussions for the empirical fitting formulae for the RM effect are presented in section 4.3. 

We analyze the RM effect of HD 17156b using our new Subaru dataset, accom- 
panied by published out-of-transit radial velocity datasets presented in F07 (Subaru) and 
Winn et al. (2009a; W09) (Keck). (The Keck dataset was originally taken by F07, and W09 
refined the radial velocities and added new data.) The out-of-transit data are needed to deter- 
mine the Keplerian orbital parameters of HD 17156b. After performing this analysis of the new 
transit data, we also repeated the analysis using all of the published radial-velocity datasets 
around the transit phase (N08, C08, and B08) for comparison. In all of our analyses we consider 
the error in each radial velocity data point to be the quadrature sum of the internal error and 
3 m s~^ (representing possible systematic error such as "stellar jitter"). 

For radial velocity fitting, we fix stellar and planetary parameters based on Barbieri et al. 
(2008) as follows; the stellar mass Mg = 1.24 [Mq], the stellar radius Rg = 1.44 [Rq], the radius 
ratio Rp/R.s = 0.0727, and the semi-major axis a = 0.1614 [AU]. We also fix the quadratic limb- 
darkening parameters for the spectroscopic band containing many Iodine absorption lines to 
Ml = 0.39 and U2 = 0.37 based on the tables of Claret (2004). In addition, we adopt an updated 
transit ephemeris based on simultaneous transit photometry as Tc = 2454777.94761 [HJD] and 
P = 21.2165298 days based on unpublished data from the Transit Light Curve project (see e.g. 
Holman et al. 2006, Winn et al. 2007). The adopted parameters are summarized in Table 2. 
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Our model has 6 free parameters describing the HD 17156 system: the radial velocity 
semiamplitude the eccentricity e, the argument of periastron cj, the orbital inclination i, 
the sky-projected stellar rotational velocity V smls-, and the sky-projected angle between the 
stellar spin axis and the planetary orbital axis A. Our model also has one free parameter, a 
velocity offset, for each independent radial velocity dataset (t^i: our Subaru dataset, V2: Subaru 
in F07, v^: Keck in W09, v^: HET in COS, v^: HJST in COS, v^: OAO in NOS, vj. TNG in 
BOS). 

We calculate the statistic 



where Wj,obs are the observed radial velocity data and ^i^caic are the values calculated based on a 
Keplerian orbit and on the RM calibration formula given above. We determine optimal orbital 
parameters by minimizing the statistic using the AMOEBA algorithm (Press et al. 1992), 
and estimate confidence levels based on the criterion Ax^ = 1.0 when a parameter is stepped 
away from its optimal value and the other parameters are re-optimized. 

4. Results and Discussions 

4.1. Results for the Key Parameter A 

We first fit our Subaru dataset with the out-of-transit Subaru (F07) and Keck (W09) 
datasets. Figure 2 shows our Subaru dataset plotted with the best-fitting model curve. We 
find A = 10.0° ±5.1°, indicating a fairly good alignment between the sky projections of the 
stellar spin axis and the planetary orbital axis. This value is similar to the result quoted by 
COS (A = 9.4° ±9.3°), which was based on F07, NOS, and COS datasets. 

Next we conduct another analysis using all of the published datasets (namely, with F07, 
NOS, COS, BOS, and W09 datasets). Figure 3 plots all radial velocity data with the best-fitting 
model curve. In this case we find A = O.S° ±4.3°, and the results for the other parameters are in 
agreement with the previous analysis. The results of both analyses are summarized in Table 3. 
We note that the central value of A for the combined analysis is somewhat biased by the TNG 
dataset. This is because the contribution to from the TNG dataset is disproportionate 
(Ax^ = 55.73 from 2S TNG radial velocities), whereas the ratio of Ax^ contribution and the 
number of data points is approximately unity for all of the other datasets. This suggests that 
the errors in the TNG dataset may have been underestimated. For this reason the results of 
the joint analysis should be interpreted with caution. 

In addition, we conduct the same fitting procedure for the OAO dataset; for the HET, 
HJST, and OAO datasets; and for the TNG dataset alone, in order to make comparisons 
with previous studies (NOS, COS, and BOS). Note that we employ the same assumptions as 
described in section 3, and we incorporate the out-of-transit velocities from the F07/Subaru 
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and W09/Keck datasets. The results are A = 57° ± 23° (as compared to A = 62° ± 25° by N08) 
for the OAO dataset; A = 11.0° ± 11.8° (A = 9.4° ±9.3° by COS) for HET, HJST, and OAO 
datasets; and A = -33.7° ±9.7° (A = -4.8° ±5.3° by BOS) for the TNG dataset. Consequently, 
we confirm that the OAO dataset indeed implies a large spin-orbit misalignment, but its poorer 
radial velocity accuracy and time-resolution would have caused the fluke with lower statistical 
significance. In addition, although our analyses are well consistent with N08 and COS, we are 
not able to reproduce the BOS result. The systematic difference may be caused by the fact that 
BOS first fitted only the out-of-transit data and determined a radial velocity offset (namely, 
the baseline of the Keplerian motion) of the TNG dataset using fewer samples than we have 
here. Such a treatment might be problematic since the RM amplitude of HD 17156b is only 
~ 15 m s~^, and a difference of offset velocity as much as several m s~^ would cause a large 
systematic difference in A. 

4-2. Assessments of Some Systematic Errors 

In the above analyses, we fixed some parameters as summarized in Table 2, which were 
independently determined by previous high accuracy photometric and spectroscopic studies. In 
order to check the robustness of our results and to estimate the level of systematic errors, we 
retry the fitting for our Subaru dataset, the F07 Subaru dataset, and the W09 Keck dataset 
with different adopted parameter values as follows; ui = 0.59 (a greater limb-darkening case); 
Ui = 0.19 (a smaller limb-darkening case); a/Rs = 21.1 (Icr lower limit in W09); and a/Rs = 25.0 
(Icr upper limit in W09). We note that we find Rp/Rg within its error reported in BOS and 
W09 does not make a significant difference in results (only V^sinJ^ is affected, and only by a 
few percent). 

We first find that /i , e, u, and offset velocities are insensitive to the fixed parameters. 
In addition, respective results for the spin-orbit alignment angle A are 10. S° ±5.6° {ui = 0.59); 
9.9° ±5.0° (mi = 0.19); 5.S°±2.9° (a/i?, = 21.1); and 13.4° ±7.0° (a/i?, = 25.0). Consequently, 
we confirm all results of A for the test cases are well within our la uncertainty. 

In contrast, we find appreciable systematic errors in the inclination i and Vsinlg. Results 
for the inclination i are S7.52°±0.36° (mi = 0.59); S6.99°±0.3r (mi = 0.19); S5.11°±0.21° 
{a/Rs = 21.1); and S7.95° ±0.40° (a/i?, = 25.0). We thus conclude that the dependence of the 
inclination on the limb-darkening parameters is small (within ~ la error), but the dependence 
on a/Rs is as large as a few degrees (over 3a). When the covariance with the a/Rg parameter 
is taken into account, our constraint on i based on the RM anomaly is approximately S4.90°- 
SS.35°. This is in good agreement with the constraint by W09 (S5.4°-SS.3°). On the other hand, 
we derive the projected velocity of stellar rotation Vsinis as, 4.0S±0.32 km s~^ (ui = 0.59); 
4.1S±0.31 km s^i (ui = 0.19); 5.0S ± 0.41 km s'^ {a/Rs = 21.1); and 4.02 ± 0.30 km s'^ 
{a/Rs = 25.0). Based on these results we find that the result for V^sinJ^ is more robust than 
that for i, but the true error is approximately 1 km s~^ (the span of the results just quoted). 
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We also consider possible time-correlated errors in the radial velocity datasets, due 
to possible instrumental or astrophysical effects. When present, this so-called "red noise" 
complicates the accurate estimation of parameter errors (Pont et al. 2006). To check for red 
noise in the Subaru dataset, we calculate the standard deviation of the residuals between the 
data and the best-fitting model, and also the time-averaged (binned) residuals. In the absence 
of time-correlated noise, the standard deviation of n-point binned residuals ((T„) would decrease 
as 



where m is the number of binned residuals. We find the standard deviation of residuals to be 
0"! = 3.42 m s~^, and if there is no correlated noise then we expect the standard deviation of 
5-point binned residuals to be (Xs cxp = 1-58 m s^^. The actual standard deviation of residuals 
with 5-point binning is cTs^act = 1-69 m s~^. This is only 7% larger than the expected value; 
apparently the level of time-correlated noise is small. 

4-3. Notes on the RM Formulae 

Winn et al. (2005) pointed out a systematic difference in the amplitude of the RM 
anomaly between the results of RM simulations and the analytic OTS formula for the HD 209458 
system. A preliminary analysis of Hirano et al. (in prep.) indicates the systematic difference 
between the OTS formula and the RM simulations is due to the fact that the former com- 
putes the spectral shift via intensity-weighted average while the latter via the cross-correlaiton. 
Specifically, they showed that the amplitude of the systematics depends on the spin rotation 
velocity of the host star. Following their procedure, we also confirmed the similar systematics 
for the HD 17156 system using the Subaru analysis routine (section 3). In addition to that, 
we show in this subsection that the input value of the spin rotation velocity, more specifically 
l^sinJs, in the empirical RM formula does not bias the estimate of the other parameters. 

When we fit our Subaru dataset using the OTS formula, we find Vsinlg = 5.03 km s~^, 
which is ~20% higher than the result using the simulation-based formula. The results for 
all the other parameters are essentially unchanged. Thus the choice of the RM calibration 
formula corresponds to a ~20% change in the inferred value of V^sinJ^,. We also expect that 
the formula of Gimenez (2006) would give similar results to the OTS formula, as both formulas 
were calculated under the same assumption: that the anomalous radial velocity is equal to the 
shift in the intensity-weighted mean wavelength of the absorption lines. 

In the previous studies, the Vsinlg value assumed for the simulation was generally in 
agreement with the resultant fitted value. However, we obtained V^sinJ^ = 4.18 ±0.31 km s~^ 
from the RM analyses even though we assumed V^sin = 2.6 km s^^ in the simulations. Ideally 
the fitted value of V sin Is should be consistent with the assumed one. In order to check the 
dependence of assumed Vsinlg value on fitting results, we reran the simulations assuming 
Vsinlg = 4.2 km s~^ for the convolution kernel instead of the original choice of 2.6 km s~^. The 
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revised empirical formula is 




) 



2' 



Av = -fvp 1.367-0.505 



V sin J, 



(4) 



s 



With this revised formula, we repeated the analysis using the Subaru dataset along with the 
F07 and W09 datasets. We find l^sinJ^ = 4.22 ± 0.31 km s~^ and the values of the other 
parameters are essentially unchanged. The results are now self-consistent and suggest that the 
choice of an empirical formula does not make a significant difference. Thus we conclude that the 
determination of A of the current system is fairly insensitive to the choice of the RM formula. 

5. Summary 

We observed a full transit of HD 17156b with the Subaru 8.2m telescope, and measured 
radial velocities with the highest time-resolution and accuracy that have yet been presented 
for transit observations of this system. We analyzed the RM effect and found that the sky- 
projected spin-orbit alignment angle of the HD 17156 system is A = 10.0° ±5.1°. This result 
supersedes the previous claim of A = 62° ±25° by N08, which was due to lower precision and 
poor statistics of the previous datasets. Thus we conclude that the HD 17156 is not strongly 
misaligned, and as such it joins the majority of other systems for which the RM effect has been 
detected (Fabrycky & Winn 2009 and references therein). If the projected spin-orbit angle is 
representative of the true spin-orbit angle, then the orbital tilt of HD 17156b is not so different 
from the orbital tilt of the Jupiter relative to the equatorial plane of the Sun. 

The other TEPs for which RM observations have been reported include HAT-P-2b and 
X0-3b. Interestingly the outcomes of the measurements of the two systems were quite different; 
HAT-P-2b is also well-aligned (Winn et al. 2007b; Loeillet et al. 2008) while X0-3b has a 
significant spin-orbit misalignment (Hebrard et al. 2008; Winn et al. 2009b). The existence of 
at least some highly-inclined planets would support the notion that some TEPs have migrated 
through planet-planet scattering or the Kozai effect and tidal friction. However, with just 
three samples one cannot draw robust conclusions (see e.g. Fabrycky & Winn 2009). Further 
observations of the RM effect for TEPs are desired. Recently the highly eccentric, long period 
planet HD 80606b has been added to the list of TEPs. Observations of the RM effect for that 
system were presented by Moutou et al. (2009), but because of their limited time coverage a 
robust determination of A is not yet possible. It will be very interesting to see the rest of the 
RM anomaly filled in with future observations. 

We acknowledge invaluable support of our observations by Akito Tajitsu, a support 
scientist for the Subaru HDS. This paper is based on data collected at Subaru Telescope, which 
is operated by the National Astronomical Observatory of Japan. The data analysis was in part 
carried out on common use data analysis computer system at the Astronomy Data Center, 
ADC, of the National Astronomical Observatory of Japan. N.N. is supported by a Japan 
Society for Promotion of Science (JSPS) Fellowship for Research (PD: 20-8141). This work was 
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the summit of Mauna Kea has always had within the indigenous Hawaiian community. 
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Fig. 1. Radial velocities taken with the Subaru HDS. The values and errors are presented in table 1. 
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Fig. 2. The upper panel: Radial velocities and the best-fit curve of HD 17156 as a function of orbital 
phase. Only Subaru and Keck datasets are used. The lower panel: Residuals from the best-fit curve. 
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Fig. 3. The upper panels: All published radial velocities and the best-fit curve of HD 17156 as a function 
of orbital phase. The upper figure shows the entire orbit and the lower figure does the zoom of transit 
phase. The lower panels: Residuals from the best-fit curve. 
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Table 1. Radial velocities obtained with the Subaru/HDS. 



lime [xlJJJJ 


Value [ni s ^\ 


Error [m s '^'\ 


z4o4 ( ( (".oUboi 


71 /I C 


01c: 
o.io 




by. do 


oc 

o.zo 


z4o4(^ ( ^ .oil (b 


oy.iy 


o.io 


z4o4r77.oi40t) 


1 A\J 


nc 


0100c 


00. z5 


o.o4 


z4o4 / ( .olibb 


Di.UO 


O.io 


z454r ( .oZbvJO 


59. DO 


O.io 


z4o4 ( ^ .ozooU 


cro no 


o.iU 


z4o4 i ( ( .C>6l\JO 


54. do 


O.io 


z454r ^o^o/^ / 


5z.5b 


o.iz 


z4o4 < (^.o4ziy 


4o.oy 


6.66 


z454r ^o4yUc> 


/I c 00 

45. oz 


O.U7 


z4o4/^ ^o5z4o 


44.00 


o.io 


z4o4 ^ /^.oooo / 


41. 94 


Q1 


O/lC/1777 Oi^ni:;^; 

z4o4777.odUdd 


T nn 


/.o7 


0/1CC/1777 cfj/inn 
z4o4 (".0041)9 


00.4/ 


Z. /O 


/i c /i 777 o^;7c;i 

z4o4/^ {.Qbibi 


c no 


n /I 
O.U4 


/i c /i 777 07nno 
z454/^ (" (.Q{[j\)6 


T 

0/.9/ 


cn 


0/1P;/1777 07/1Q!^ 

z4o4 ( ( .Qi 4^0 


Q 1 /in 


Q Of? 
O.ZO 


0/iC/i777 07770 


00.16 


000 

z.oo 


0/1F;/1777 oQion 
z4o4 ( ( ( .CtoiZD 


QQ QO 
00. OZ 


o.Uo 


/i C /i 777 OO/li^l 

z454/^ (" ^oo4Di 


00 1 c 

oz.io 


00 


/i c /i 777 ooono 
z454/^ (" ^oooL)^) 


0(3. 6z 


1 /I 
o.i4 


0/iC/i777 oni /1C 

z4o4r77. 59145 


4o.oU 


2.74 


0/iC/i777 On/107 

z454r7r.oy4o/ 


00 00 
00. OZ 


00 


0/ic;/i777 cncon 
z4o4 ( ( .oyozy 


on 00 

oa.zz 


Q K 

z.oo 


/i c /i 777 nni 71 
z454r ( (^.yiJi / i 


on /in 
09. 4U 


1 n 
O.iU 


0/ic:/i777 nnci a 
z454/^ (" ^yL)5i4 


01 1/1 
0I.I4 


000 
2.80 


0/iC/i777 nnoc^; 
z4o4r77.yL)ooD 


on 00 
zy.zz 


on 
0.2U 


0/iC/i777 ni ino 

z454r r ^yiiyo 


01 no 

ol.9o 


on 

o.oU 


0/1 777 01 !^/1 n 


zy.D4 


0.00 


2454777.91882 


25.67 


3.11 


2454777.92224 


24.40 


2.81 


2454777.92566 


19.20 


3.32 


2454777.92909 


17.00 


3.01 



11 



Table 1. (Continued.) 
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Q 1 K 

o.io 


z4o4/^ (^.yDDO/ 


-zU.o4 


O OO 

6.6Z 


z4o4 ^ ^ .y / uzy 


QQ OC 

-oo.zo 


Q y1 S3 


Z4:b4:{7 I .\) i6 i I 


OO 1 n 


O 'TO 


z454r r^y M 14 


oo no 


o 1 n 


z4o4 (".yoUo / 


-0( .OKJ 


Q 70 
O. /O 


z4o477/.yooyy 
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o 1 n 
O.iy 


z4o4 / / .yo 4d 


/1 nc 


Q 7n 


z4o4r77.yyL)o5 


/I /I TO 

-44. 


O.DY 


z4o4 / 77.yy4zo 


o c 1 

-oo.ib 


cS.71 


z4o4 (".yy / 


A^ dQ 


o.4i 


z4o4 /o.UUiiz 


-4b. zo 


o 1 n 


OAKATTQ r\r\AKA 


on OA 
-oy.z4 


O K 

o.oo 


z4o4 r /o.UU /yb 
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o /I c /I TTo n 1 1 on 
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oa nn 
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00 A'7 

-66 Ai 


Q on 
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o oo 
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z4o4/^ /^o.Uzibo 


/I O 1 c 


o.4i 


z4o4 i ^ o.UzoUo 




Q 01 


z4o4r7o.L)zooi 


AO T^Cl 


O 'TO 


z4o4 ('o.Uoiy4 


A'7 07 

.6i 


Q Oyl 


z4o4/^ /^O.UJo/L) 


-4b. ob 


o o c 

o.oo 


z4o4/^ /^O.Ujyo^) 


CI oo 


o oo 


z404 ( ^ o.U4oyD 


-00. Uo 


0.10 


2454778.04809 


-54.08 


3.04 


2454778.05222 


-57.89 


3.31 


2454778.05635 


-57.39 


3.01 


2454778.06048 


-55.33 


3.39 
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Table 1. (Continued.) 



z4o4 ( /o.UD4Di 


-oy.zo 


O.Zi5 


9/1 778 nfiS7/l 


-Do.oO 


Z.ol 


2454778.07288 


-69.28 


3.00 


2454778.07701 


-74.63 


3.42 


2454778.08114 


-72.57 


2.90 


2454778.08527 


-80.46 


3.12 



Table 2. Adopted stellar and planetary parameters. 



Parameter 


Value 


Source 


Ms [Mq] 


1.24 


BOS 


Rs [Rq] 


1.44 


BOS 


Rp/Rs 


0.0727 


BOS 


a [AU] 


0.1614 


BOS 


Ul 


0.39 


Claret (2004) 


U2 


0.37 


Claret (2004) 


Tc [HJD] 


2454777.94761 


Winn et al. in prep. 


P [days] 


21.2165298 


Winn et al. in prep. 
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Table 3. Fitted values and uncertainties of the free parameters. 





Subaru and Keck 


all datasets 


Parameter 


Value 


T T J. * J. 

Uncertainty 


Value 


Uncertainty 


A [m s J 


275.99 


±1.27 


277.77 


±1.19 


e 


0.6801 


±0.0019 


0.6835 


±0.0017 


roi 


121.62 


±0.42 


121.73 


±0.41 


■ rol a 

I [J 


87.21 


±0.31 


86.95 


±0.27 


V siiiis [km s J 


4.18 


±0.31 


4.07 


±0.28 


\ rol 

A [°] 


10.0 


±5.1 


0.8 


±4.3 


r —11 
vi [m s J 


97.46 


±1.57 


98.23 


±1.51 


r , -11 




4-0 1 n 


1 
ao. iz 


4-0 no 


W3 [m s "'"] 


-0.11 


±0.64 


-0.20 


±0.64 


114 [m s^"'"] 






74.21 


±2.71 


Us [m s""'"] 






-7850.55 


±3.01 


vq [m s""'"] 






142.34 


±3.29 


v-j [m s""'^] 






146.68 


±1.89 



v\ ^ V7 are offset velocites for respective datasets (see Section 3). 
^: Systematic errors are not included in the uncertainties (see Section 4.2). 



14 



References 



Barbieri, M., et al. 2007, A&A, 476, L13 
Barbieri, M., et al. 2008, ArXiv e-prints (BOS) 
Barker, A. J. k Ogilvie, G. I. 2009, ArXiv e-prints 

Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C., Dosanjh, P., & Vogt, S. S. 1996, PASP, 
108, 500 

Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580 
Claret, A. 2004, A&A, 428, 1001 

Cochran, W. D., Redfield, S., Endl, M., & Cochran, A. L. 2008, ApJL, 683, L59 (C08) 

Fabrycky, D. C. k Tremaine, S. 2007, ApJ, 669, 1298 

Fabrycky, D. C. k Winn, J. N. 2009, ArXiv e-prints 

Fischer, D. A., et al. 2007, ApJ, 669, 1336 (F07) 

Gaudi, B. S. k Winn, J. N. 2007, ApJ, 655, 550 

Gimenez, A. 2006, ApJ, 650, 408 

Hebrard, G., et al. 2008, A&A, 488, 763 

Holman, M. J., et al. 2006, ApJ, 652, 1715 

Ida, S. & Lin, D. N. C. 2004, ApJ, 616, 567 

Kurucz, R. L., Furenlid, I., Brault, J., & Testerman, L. 1984, in Solar flux atlas from 296 to 1300 nm, 

ed. R. L. Kurucz, I. Furenlid, J. Brault, & L. Testerman 
Lin, D. N. C, Bodenheimer, P., & Richardson, D. C. 1996, Nature, 380, 606 

Lin, D. N. C. & Papaloizou, J. 1985, in Protostars and Planets II, ed. D. C. Black & M. S. Matthews, 
981-1072 

Loeillet, B., et al. 2008, A&A, 481, 529 

Marzari, F. & Weidenschilhng, S. J. 2002, Icarus, 156, 570 

McLaughlin, D. B. 1924, ApJ, 60, 22 

Moutou, C, et al. 2009, ArXiv e-prints 

Nagasawa, M., Ida, S., & Bessho, T. 2008, ApJ, 678, 498 

Narita, N., et al. 2007, PAS J, 59, 763 

Narita, N., Sato, B., Ohshima, O., & Winn, J. N. 2008, PASJ, 60, L1+ (N08) 

Noguchi, K., et al. 2002, PASJ, 54, 855 

Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118 

Pont, F., Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C. 

The art of scientific computing (Cambridge: University Press, — cl992, 2nd ed.) 
Rasio, F. A. & Ford, E. B. 1996, Science, 274, 954 
Rossiter, R. A. 1924, ApJ, 60, 15 

Sato, B., Kambe, E., Takeda, Y., Izumiura, H., & Ando, H. 2002, PASJ, 54, 873 
Takeda, G. & Rasio, F. A. 2005, ApJ, 627, 1001 
Winn, J. N., et al. 2005, ApJ, 631, 1215 

Winn, J. N., Holman, M. J., & Fuentes, C. I. 2007a, AJ, 133, 11 
Winn, J. N., et al. 2007b, ApJL, 665, L167 

15 



Winn, J. N., et al. 2009a, ApJ, 693, 794 (W09) 

Winn, J. N., et al. 2009b, ArXiv e-prints 

Wu, Y. k Murray, N. W. 2003, ApJ, 589, 605 

Wu, Y., Murray, N. W., & Ramsahai, J. M. 2007, ApJ, 670, 820 



16 



